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Abstract. By varying the orientation of the applied magnetic field with respect to the normal of a two- 
dimensional electron gas, the chemical potential and the specific heat reveal persistent spin splitting in all 
field ranges. The corresponding shape of the thermodynamic quantities distinguishes whether the Rashba 
spin-orbit interaction (RSOI) , the Zeeman term or both dominate the splitting. The interplay of the tilting 
of the magnetic field and RSOI resulted to an amplified splitting even in weak fields. The effects of changing 
the RSOI strength and the Landau level broadening are also investigated. 

PACS. 71.70.Di Landau levels - 71.70.Ej Spin-orbit couphng - 73.20.At Electron density of states 



1 Introduction 

New physics has been exhibited by two-dimensional elec- 
tron gas (2DEG) systems even at weak magnetic fields. 
Transport properties become highly anisotropic with the 
direction of the applied magnetic field B playing a cru- 
cial role [1112,3,4,5,6 . The anisotropy is considered to be 
a product of the interplay of the spin-orbit coupling, the 
Zeeman splitting and the orientation of B with respect to 
the 2DEG plane among others. These factors determine 
the density of states (DOS) of the system. As many ob- 
servables, such as the conductance are proportional to the 
DOS at the Fermi energy [7), tuning each of these char- 
acteristics is indispensable in the design of devices that 
require the manipulation of spins. Such devices are being 
eyed in spintronics and quantum information systems j8j. 

The spin-orbit interaction experienced by the electrons 
in a 2DEG can be due to the structure inversion asymmet- 
ric potential (Rashba) or the bulk inversion asymmetry 
(Dresselhaus) . Both lead to a zero- field spin splitting. The 
former is of particular interest due to its versatility as it 
can be controlled by the gate voltage pilTU]. Moreover, in 
some III-V [n], 11- VI [TO] and Si-based [12] semiconduc- 
tors the Rashba spin-orbit interaction (RSOI) normally 
dominates. Beating patterns observed in transport prop- 
erties relative to the applied magnetic field were attributed 
to the RSOI [HldnillS. 

Tilting the magnetic field with respect to the 2DEG's 
normal axis has revealed a rich variety of physical signa- 
tures. For example, it recovers the resonant spin Hall ef- 
fect after being suppressed by an impurity scattering 0]. 
Depending on the angle, a field's orientation can destroy 
some fractional quantized Hall states or turn anisotropic 
phases into isotropic ones [S]. It also causes the ringlike 
structures in the longitudinal resistivity versus magnetic 



field to collapse ^. With such interesting tilt- induced ef- 
fects, this work focuses on the influence of the applied 
magnetic field's direction on the thermodynamic proper- 
ties of a 2DEG with Rashba and Zeeman interactions. 
We now investigate how a finite in-plane component of B 
modifies the chemical potential and the specific heat. The 
shape of their oscillations will be shown to depend on both 
the magnitude and orientation of B. Moreover, robust 
spin splitting is observed even at the weak field region. 
The behavior of the thermodynamic quantities will be ex- 
plained through the competing dominance of the Rashba 
and Zeeman interactions in various field regimes. 



2 Eigenvalues with RSOI and tilted magnetic 
fields 

The 2DEG considered here is lying on an x — y plane 
and is subject to an external magnetic field of magnitude 

b^Jbi + bi 



IBI +BI + B^. The tilt is the angle 9 that B makes 

with the z-axis. The Hamiltonian of a single electron in 
this system can be written as 
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where the first term is the free particle energy, the second 
is the Rashba spin-orbit interaction, and the third is the 
Zeeman energy. The strength of the RSOI, assumed to be 
constant for this present work, is indicated by the param- 
eter a. Here m* is the electron's effective mass, cr are the 
Pauli matrices and k is the wave vector. The magnitude 



of the latter is determined by A; = 



'rA, where e is 



the electronic charge, h is Planck's constant over 27r, and 
A is the magnetic vector potential. The magnetic moment 
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is /X = ^gfiBcr where g is the Lande factor and is the 
Bohr magneton. Although the g-factor varies with B T4', 
[T5] . it was shown that the S-dependence of the 2DEG 
specific heat is not significantly altered [U]. Our focus 
now is on the effects of the RSOI and tilted fields. For 
this purpose, we keep g constant, that is, g — 2 for an 
electron. 

In solving the Schrodinger equation, 

/ Hn + E iTa + \ / Y.n=0 On^n \ - () 

(2) 

we express the wave function solutions in terms of a su- 
perposition of the harmonic oscillator basis functions 0„ 
where the coefficients a,i and 6„ are the unknown spin- 
up and spin-down complex coefficients of the nth Lan- 
dau level, respectively. Here iJ„ = \{aa^ -\- a^a), Qj = 
^BBj/huJc is the j-th component of the rescaled Zeeman 
energy, T — 2a^ C/h^c is the rescaled Rashba parame- 
ter, = m* /2h? and fl± = fix iifly In this work, energy 
units are given in terms of the cyclotron energy HlJc where 
ojc — eBz/m* . Equation ^ yields a set of secular equa- 
tions 

(^n + ^ + [2^ - a„ + iYy/n + 1 6„+ 1 + 6„ = 0, (3) 



- iT^/n ttn-i + n+Un + yn + - - - E j bn = 0. (4) 

Note that the secular equations derived by Rashba [TB] 
are recovered when B is normal to the 2DEG plane, that 
is, when we let 12+ = i7_ = 0. 

Taking advantage of the Landau level crossing, a solu- 
tion can be obtained by invoking that opposite spin states 
of adjacent Landau levels En are equally probable. The 
eigenvalues Ei^ can be shown to be 

Et =n+^ + ^{V^^+V^i) 

±iy^[r(^Am- V^) + 2r2j%4i2+r2_. (5) 

This result will be used to determine the density of states 
(DOS) which, in turn, will be used in calculating for other 
thermodynamic quantities. 

An important caveat has to be issued regarding ([5]) . It 
should not be evaluated for the perpendicular case 9 = 0. 
For such a case, J7± has to be set to zero right at the 
secular equations (|3]) and (Hj). 

The eigenvalues of equation ([5]) do not exhibit cross- 
ings with respect to increasing RSOI strength. This is un- 
like in the case when B is perpendicular to the 2DEG 
plane. The absence of intersections was conclusively at- 
tributed to the in-plane component of B [T7]. How these 
findings subsequently affect the chemical potential and the 
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Fig. 1. The chemical potential as a function of the magnetic 
field. Here a = L5 x 10"" eV-m, T = 0.5 meV and 6 = 
30°. We used T = 10 mK which is of the same order as in 
experiments [2 • The integers coinciding with the dips indicate 
filling factors. The Roman numerals indicate the different B 
regimes as discussed in section [S] 
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Fig. 2. Comparison of the different energy contributions. Here 
61 = 30°, ai = 1.5x10"" eV-manda2 = 5.5x10"" eV-m. The 
Roman numerals indicate the different B regimes as discussed 
in section |3l 

specific heat of a 2DEG will be analyzed in this paper. It 
will be shown here that spin splitting in the weak B region 
can be greatly smeared in the DOS but it can be clearly 
shown by the chemical potential and the specific heat. 



3 Chemical Potential 

The chemical potential /x is the change in energy of a ther- 
modynamic system if an additional particle is introduced, 
with the entropy and volume held fixed |18| . At zero tem- 
perature T, the chemical potential fi is the Fermi energy. 
When T is nonzero, /i varies with T and its behavior is 
shown in reference [TO] for a 2DEG in a perpendicular 
magnetic field. 
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Fig. 3. The chemical potential for different tilt angles (top 
frame) and their corresponding DOS for E = 15.8 meV (bot- 
tom frame). Here T = 10 mK, F — 0.5 meV and a = 
1.5 X IQ-" eV-m. 
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Fig. 4. The chemical potential for different Rashba interaction 
strengths (top frame) and their corresponding DOS for E — 
15.8 meV (middle and bottom frame). Here T = 4.2 K, T = 0.5 
meV and 9 — 30°. The dotted line corresponds to 02 ~ 5.5 x 
10"^^ eV-m while the solid line to qi = 1.5 x 10"" eV-m. 



Experiments on 2DEG devices are usually constrained 
to a constant particle density [31 IT51I20] . Hence, to simulate 
conditions similar to experiments, the electron concentra- 
tion is set to = 3.2 X lO^^cm"^. This is of the same order 
of N in other research works [2l[51[51 [Tn[ni[n] . Knowing N, 
fi can be numerically derived from 



N 



f{E)BOS{E)dE, 



(6) 



where f{E) = l/(exp[(£'- /x)/fcBT] + 1), ks is the Boltz- 
mann constant and DOS(i?) is the density of states. Since 
most experiments reveal a broadenend DOS fE,10,14,lS] 
and theoretical fits [TSII^ agree well with a Gaussian 
form, we also utilize it here. The DOS can be expressed 



as 



BOS{E)^'-^Y. 



1/2 



— exp 



{E-E, 



±\2 



2r2 



(7) 

where 7^ is the broadening parameter. The values F = 0.1 
meV and 0.5 meV are used in this work which are typical 
for numerical simulations in 2DEG [21] • In evaluating E^ 
in ([5]), we used to* = 0.05me where rrie is the free electron 
mass. This is the effective mass of a 2DEG in InGaAs 
heterostructures [9 . For simplicity we let By — 0, = 
B sin 9 and B^ = B cos 0. Equation ^ is solved with a 
maximum percent error of 10~^% in N. 

In figure [1] we study the chemical potential oscilla- 
tions as a function of B. The oscillations as B is increased 
are due to the depopulation of the energy levels. A good 
measure of the occupation of states is the filling factor 
v. It is the nominal number of filled energy levels and is 
inversely proportional to the component of the magnetic 
field normal to the 2DEG plane, that is, v — hN/eBz- 
In the quantum Hall effect, the conductivity is quantized 
when v takes on integer values. Because the prefactor of 
the DOS in (O is proportional to B^, when B is increased, 
for a fixed 9, the degeneracy of the energy levels also in- 
creases. Electrons from the last occupied level will then 



move to lower levels until that former level becomes com- 
pletely depopulated. This depopulation is marked by a 
drop in fi. This is clearly illustrated in figure [1] by the lo- 
cations of the dips that fall right at integer values of the 
filling factor v. 

It can be observed that /i behaves in three distinct 
ways according to its shape versus B as shown in figure 
[T] The description will be presented here and the next 
paragraphs will be devoted to its discussion. At weak B 
{v > 50), the oscillations resemble spikes of approximately 
uniform widths with the even and odd integer fillings hav- 
ing similar heights (top frame). As B increases, the shape 
of /i gradually changes. When B is moderately strong 
(8 < < 50), alternating spikes appear where one is 
broader and taller and centered at even integer v while the 
other is thinner and shorter and centered at odd integer 
v (middle frame). As B further increases, fi again slowly 
changes shape until finally, at very strong B [v < 8), 
large triangular waves each with a steep drop at odd v 
develop. The triangular waves each has a short but sharp 
spike that is centered at even v (bottom frame). The three 
regions distinguished here are, respectively, identified as 
follows: (I) the regime where Rashba SOI dominates, (II) 
the region where the Rashba and Zeeman interactions are 
comparable, and (HI) where the Zeeman term prevails. 
Ascribing region I (HI) as RSOI (Zeeman) dominated is 
based on the fact that the RSOI (Zeeman energy) is inde- 
pendent (linearly dependent) on B. These divisions based 
on the strength of the magnetic field are observed in all fi 
data obtained in this work. 

When only the cyclotron energy Ec is considered, that 
is, spin is neglected, fx peaks are large only for even v 
f2?ll^ . The presence of odd v peaks for all B region in fig- 
ure [U for example, is evidence of a substantial spin split- 
ting. Although the Zeeman energy Ez also contributes to 
the appearance of odd v peaks, it is negligible in the weak 
B region (Region I). To track the contributions of the two 
sources of splitting, we compare i?c, Ez and the RSOI en- 
ergy contribution Eb,- Refer to figured] We considered two 
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values of the RSOI parameter a tliat are seen in literature. 
We first tackle the case when the RSOI strength is equal to 
ai. In the region indicated by I, because Ec > Ef> > Ez, 
the splitting is mainly contributed by RSOI. In region II 
the splitting is a combined effect of Er and Ez since the 
two are comparable in this J5-range. Here /x behaves in a 
fashion where regions I and III shapes blend. In Region III 
Ec > Ez > En ■ Thus, the big peak of ^ at odd v in region 
III is dominated by Ez- The latter observation is similar 
to the integral quantum Hall effect where the plateau of 
the Hall conductivity Uxy at even v is determined by the 
large E,. while the corresponding plateau for the odd v is 
more influenced by Ez [23] . When stronger RSOI is con- 
sidered, a2 for example, the spin splitting is completely 
dominated by RSOI and Ez can be neglected. 

The distinguishing mark of figure [T] in contrast with 
previous results [5T] is the appearance of a well-defined 
spin splitting in the weak B case (Region I). This can 
be attributed to the RSOI plus the tilting of B. The 
RSOI with perpendicular B is normally associated with 
the oscillations in the magnetization pi[T5] and longitudi- 
nal resistivity [51ll0|. They have beats whose amplitudes 
are small at weak B relative to their values at high B. 
Figure [T] suggests that tilting B amplifies the oscillations 
as B ^ 0. 

The sharpness of the oscillations are unexpected since 
we made use of a Gaussian DOS. Refer to the lower frames 
of figures [3l HI and [5] which are evaluated at E « 15.8 
meV. This value is chosen since it lies below but close to 
the Fermi energy where the occupancy of f{E) sa 1. For 
E below this value, we obtain similar qualitative behav- 
ior. Closely inspecting DOS, we see the close correlation 
of jjL and DOS. We find that although the Gaussian func- 
tion remains the same, the interplay of the RSOI, weak 
B and tilting effectively changed the bell shape. In this 
case, the contributions from different E^ in the summa- 
tion in ([7]) are not uniform and at certain values of -B, all 
contributions drop. 

In figure [31 the effect of varying the direction of B 
is displayed. The data shown are limited to the weak B 
region since all /i data obtained are of the same qualitative 
behavior as figure [T] The angle 9 changes the DOS in two 
ways. First is through E^ which determines the location 
of the DOS peaks. See equation ([7|). The energy levels Ef^ 
in ([SJ , in turn, are affected by 9 through the components of 
the rescaled Zeeman energy Qj . We can conclude here that 
9 changes the frequency and phase of the ji and DOS via 
the latter's dependence on the energy levels. Second, the 
DOS has a prefactor that is proportional to Bz = B cos 9. 
The direction of B then affects the amplitude of these 
thermodynamic quantities because of the DOS prefactor. 

Figure[3|shows that increasing the RSOI strength raises 
/X. Here T — 4.2 K is chosen so that we can resolve better 
the difference between the behavior of /i for the two differ- 
ent RSOI strengths. From its definition, a larger /i signifies 
that a system requires higher energy for an addition of a 
particle to be possible. We can see from the figure that 
more asymmetric quantum wells demand larger amount 
of energy when a change in electron number is desired. 
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Fig. 5. The chemical potential for different broadening param- 
eters (top frame) and their corresponding DOS for E ~ 15.8 
meV (middle and bottom frames). To observe clearly the solid 
line, we magnify it in the bottom frame. Here T = 10 mK, 
B = 30° and a = 1.5 x 10"" eV-m. 



Less a shifts both DOS and /i to lower values. However, 
unlike the DOS, persistent spin splitting is manifested by 
[i{B) data as B ^ (cf. dotted line in figure[l|). The DOS 
spin splitting is more robust when a is larger but the cor- 
responding amplitudes of the oscillations are shorter. Here 
we can see a trade-off of broadening the Landau levels in 
exchange to a resolved spin splitting. 

The \i peaks in all our simulations are sharp unlike in 
previous results where the spin effects were ignored [35]. 
The oscillating behavior of the 2DEG with respect to B 
is better resolved in /i than in the DOS as i? — ?► 0. 

Furthermore, we show the effect of increasing the Lan- 
dau level broadening in figure [5] An increase in the broad- 
ening parameter changes the distribution of occupied states 
over the energy. Formerly unoccupied states become oc- 
cupied. This is at the expense of the once highly occupied 
states. Hence, the total number of states remains con- 
served. A large F indicates a system with more disorder 
jl9) . The latter can take the form of structural disorder 
or impurity scattering in which charged impurities inter- 
act with an electron via the screened Coulomb interaction 
[IS] . Ordinarily F does not drastically alter the qualitative 
behavior of the DOS. For a fixed A^, a small E is associated 
with a narrow DOS curve having a tall amplitude. This is 
in contrast to the case of a wide E where the DOS curve 
is broadened but its amplitude is relatively shorter. How- 
ever, based on equation (I7|), the DOS has a complicated 
dependence on both E and B, where the latter appears 
nontrivially in En through equation ([S]). Hence, the DOS 
shape in figure [5] is a result of the combined effects of 

the RSOI and the tilted field. It is therefore difficult 
to attribute conclusively to a single parameter a specific 
behavior of the DOS. 

Despite the not-so-simple dependence of the DOS on 
the different parameters, figure [5] suggests that a more dis- 
ordered system is more resistant to changes in the number 
of particles in the 2DEG resulting to a higher /i. Here we 
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find that even at a weak a small F is sufficient to dis- 
play the fi behavior expected only at strong fields. This 
is depicted by the shape of /i for F — 0.1 meV in the 
top frame of the figure [5] which is similar to region III of 
figure [H Comparing figures S] and [SJ we find that increas- 
ing the RSOI strength and the Landau level broadening 
have a similar effect on /i. Increasing the RSOI strength 
increases spin splitting which causes adjacent Landau lev- 
els to overlap. Similarly, two adjacent Landau levels will 
overlap if F is sufficiently wide. 



4 Specific Heat 

The specific heat capacity Cv is defined as the amount 
of energy required to raise a system's temperature by a 
degree while the volume is kept fixed [TB]. It is important 
to determine Cv so that engineers will know the heat tol- 
erance of devices such as those intended to be resistant to 
changes in T. 

At constant volume, Cy is given as 

= ^ 1 f{E) (E - ^) DOS(i?) dE. (8) 

The partial derivative in ([SJ will yield two terms: one 
from df/dT and another from dfi/dT. The contribution 
of these two terms are taken into account in figure [6] since 
varies with both B and T. 

Experimental and numerical data show that a 2DEG 
under a perpendicular B has a Cy that oscillates as a 
function of B with the dips coinciding with integer v [131 
[TSU^Dll^l^ . These minima coinciding with integer v are 
also demonstrated in figure [6] 

Similar with the categories of behavior depending 
on the dominant interaction, Cv{B) also has three distinct 
features. In weak fields (region I) , the peaks are composed 
of spin split levels of heights slowly increasing (see inset 
of figure El) . When B is moderately strong (region II) the 
oscillations are composed of peaks of alternating heights. 
Finally, when B is strong, wide and tall peaks emerge with 
a sharp dip that cuts through each peak. The shape of 
the Cv{B) at region III resembles that of the spinless case 
[22] . The Zeeman term in this region, creates an additional 
dip in the previously spin-degenerate Landau levels. The 
latter finding was also observed in earlier works |24j where 
the additional dip was attributed to Coulomb interaction 
in the presence of a disorder potential. 

The prominence of the odd ly at weak and moderately 
strong B gives evidence to the influence of the RSOI at 
these B regions. The same argument used in the preced- 
ing section applies here. Usually only the even v are pro- 
nounced with small dips at odd v. This common trend 
shows that the cyclotron energy Ec is usually larger than 
any spin effect [53] . The comparable size of the peaks at 
strong B implies that at this region, the Zeeman splitting 
is also comparable to the cylotron energy (cf. figure [2]). 
Both fi in figure [1] and Cy in figure [5] manifest an ampli- 
fied spin effect in all B spectrum. The plots also highlight 




beside the dips indicate the filling factors. The inset is for the 
weak B case. 

the area where one of the two spin interactions dominate 
or where both become comparable. 

The temperature used in the simulation of figure |6] is 
T = 4.2 K, the same T used in the Cy measurements 
of Ref. . We chose this temperature because the mil- 
likelvin range has been found to freeze disorder effects and 
it constrains the system to behave as an ideal free electron 
gas [22 . This was also observed in reference [55] although 
their consideration was without RSOI. The behavior of 
Cv with T is the subject of another paper. 

Another notable observation is the nonmonotonic in- 
crease in the height of the Cy peaks as B increases where a 
distinct turning point can be found around B « 1.2 T. See 
arrow pointing to — 13 in figure [G] Moreover, the dip at 
this point is displaced upward compared to Cy ~ for the 
rest of the dips. This could be reminiscent of beat nodes as 
observed in reference |13j that appears to be damped here. 
Although not shown here, we observe that the dip shifts 
to higher B with 9. This is similar to the observation that 
tilting the applied field shifts the beats of the i?-dependent 
magnetization [5] and longitudinal resistivity [3]. 

5 Conclusions 

In this work, we have shown the distinct behavior of the 
chemical potential /i and the specific heat capacity Cy of 
a two-dimensional electron gas subject to a tilted mag- 
netic field B with varying degrees of energy-level splitting 
interactions. These thermodynamic properties were ob- 
tained assuming a Gaussian density of states DOS whose 
energy levels are obtained analytically in the presence of 
the Rashba spin-orbit interaction RSOI and the Zeeman 
effect. The RSOI dominates at the weak B region while 
the Zeeman energy rules at the strong B part. Interesting 
fj, and Cy behavior occur when these two contributions 
are comparable to each other. The combined effect of the 
RSOI, the Zeeman splitting and the tilting of B yielded 
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here a damping of beating patterns in the oscillations of 
Cv and the DOS as observed in other studies found in 
literature. 

Among the tunable parameters investigated, only the 
tilt angle 9 effectively changed the phase and frequency 
of the /I oscillations. This outcome can be traced to the 
dependence of the DOS on 9 . Finally, the interplay of 
the RSOI and the tilted magnetic field resulted to an am- 
plified spin splitting as i? — >■ which was not observed 
before. The persistent splitting is demonstrated by the 
pronounced /i and Cy peaks at odd filling factors which 
is not observed in the spin degenerate case. This signifies 
that in addition to tuning the gate voltage, in varying the 
RSOI strength, the orientation of the field relative to the 
electron gas plane provides control to the magnitude of 
the spin splitting which is important in the spin filtering 
design of 2DEG devices. 
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